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ABSTRACT 

We have observed the Class I protostar TMC-lA with Atacama Millime¬ 
ter/submillimeter Array (ALMA) in the emissions of and (J = 2 — 1), 
and 1.3-mm dust continuum. Continuum emission with a deconvolve size of 
Cf.'SO X (y.'37, perpendicular to the ^^CO outflow, is detected. It most likely traces 
a circumstellar disk around TMC-lA, as previously reported. In contrast, the 
C^®0 a more extended structure is detected in C^®0 although it is still elongated 
with a deconvolved size of 3'.'3 x 2'.'2, indicating that C^®0 traces mainly a flat¬ 
tened envelope surrounding the disk and the central protostar. C^^O shows a 
clear velocity gradient perpendicular to the outflow at higher velocities, indica¬ 
tive of rotation, while an additional velocity gradient along the outflow is found 
at lower velocities. The radial prohle of the rotational velocity is analyzed in 
detail, hnding that it is given as a power-law oc with an index of ~ 0.5 at 
higher velocities. This indicates that the rotation at higher velocities can be 
explained as Keplerian rotation orbiting a protostar with a dynamical mass of 
0.68 Mq (inclination corrected). The additional velocity gradient of C^®0 along 
the outflow is considered to be mainly infall motions in the envelope. Position- 
Velocity diagrams made from models consisting of an infalling envelope and a 
Keplerian disk are compared with the observations, revealing that the observed 
infall velocity is ~0.3 times smaller than free fall velocity yielded by the dynami¬ 
cal mass of the protostar. Magnetic helds could be responsible for the slow infall 
velocity. A possible scenario of Keplerian disk formation is discussed. 

Subject headings: stars: circumstellar matter — stars: individual (TMC-lA) — stars: 
low-mass — stars: protostars 
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1. Introduction 

Protoplanetary disks play key roles in the evolution of young stellar objects (YSOs). 
In the early stages of star formation, mass and angular momentum are transferred to the 
central stars through protopla netary disks. In the later stages, these disks can become the 


sites of planet formation (e.g. 


Howard 


20131) . Protoplanetary disks are ubiquitous around 


classical T-Tauri stars (CTTSs) or Class n sources and have been observed in thermal dust 


the 1990s. Keplerian rotation is revealed for manv of the disks ( 

Guil 

oteau & Dutrev 

1998: 

Guil 

oteau et ah 1999: Simon et ah 2000; Oi et ah 

2004; Hughes et ah 

2009 


Rodriguez et al. 

2010 

). In the evolutionary stages earlier than T-Tauri stars, accreting protostars (Class I or 


0 sources) are associated with disk-like structures that exhibit continuum and line emissions. 
Thes e disk-like structures often show infall motions with little rota t ion (e .g. 


2009 


Enoch et ah 


2009 


2011 


Brinch et ah 


2009 


Yen et ah 


2010 


Jprgensen et ah 


20111) and are called 


‘pseudo disks” (Galli and Shu, 1993ab), although pseudo disks have not been hrmly 


have also revea 

ed Keplerian disks around Class 0 anc 

Class I protostars 

(Takaknwa et al. 

2012; 

Yen et al. 

20ls3; 

Harsono et al. 

2014; 

Choi et al. 

2010; 

Tobin et al. 

2012b 

), that are 


still deeply embedded in protostellar envelopes. 

Despite the ubiquity of Keplerian disks in various protostellar stages, their 
formation process from protostellar envelopes is still controversial both theoretically 
and observationally. Theoretical calculations of collapsing isother mal spheres with ide al, 


non-magiietized conditions, such as “inside-out collapse” models (iTerebey et ah 


Shuetah 


19841 : 


19871) . provide a conventional picture of star and disk formation. According to 


these models, a central protostar forms with a Keplerian disk at the center of a slowly 
rotating dense molecular cloud core as the gas begins dynamical infall, which conserves its 
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angular momentum. Larger angular momentum is gradually transferred from the outer 
region of the dense core and is re-distributed in the disk as time goes on. There are some 
observational studies conhrming this tendency that a Keplerian disk develops later so that 


the rotatio nal velocity and siz e of the dis 


YSOs (e.g. 


Ohashi et al. 


1997 


Yen et al 


r bec ome larger in later evolutionary stages of 


20131 ). Although these studies observationally 


suggest the growth of Keplerian disks in the course of YSO evolution, the dynamics of 
each source are not spatially resolved in these studies. This prevent us from directly 
investigating how a disk develops in a protostellar system. Recent theoretical simulations 
for disk formation around protostars, including MHD calculations, are still under heated 


discussion. It is reported that such simulations have difficulties in trans forming enve 
to Ke plerian dis ks because of s t rong r nagne t ic braking or fr a gmen tation flMellon fc Li 


opes 


2008 


20091 ). However, 


Tomida et al. 


(12013 


20151) : 


Machida et al. 


fl201l[) point out that non-ideal 


MHD effects can suppress angular momentum transport by magnetic helds a nd enable 


form ation of circumstellar disks in the early phase of star formation (see also 


Joos et al. 


20121 ). The formation of Keplerian disks is closely related to the dynamical transition 


from the infall motion in protostellar envelopes to the Keplerian rotation. In fact, such a 


transition is being observationallv revea 

ed around several 

rrotostars bv the most recent 

studies (’ 

Vlurillo et al. 

2013; 

Chou et al. 

2014; 

Ohashi et al. 

2014; 

Takakuwa et al. 

2014; 


Yen et al. 


20141) . The number of samples is, however, still too limited to perform any 


statistical study. In order for us to unveil the disk formation mechanism, it is essential to 
observe the two kinds of dynamics (infall and rotation) around protostars and increase the 
number of direct detections of the transition between the two. 


In order to investigate the expected transition from the infall motion to Keplerian 
rotation, we made observations of the Class I protostar TMC-lA (IRAS 04365-1-2535) with 
the Atacama Large Millimeter/Submillimeter Array (ALMA) in J = 2 — 1 (230.5380 
GHz), J = 2 — 1 (219.5604 GHz), and 1.3-mm continuum emission. As its name 
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suggests, our target TMC-IA is in the Taurus molecular cloud, which is one of the closest 
star-forming regions {d = 1 40 pc). Its bolometric luminosity -f^boi — 2.7 Lq and bolometric 


temperature Tboi = 118 K flKristensen et al 


20121) indicate that TMC-lA is typical Class 
I protostar. The sys temic velocity of TMC-lA was measured at 14sr = 6.4 km s“^ by 


Terebev et al. 


(119891 ) 


7 

^ molecular on 

thow (tdvr> ~ 2.5 X 10 

^ vr) 

1995; 

Tamura et al. 

1996; 

Chandler et al. 

1996) 


yr) is associated with TMC-lA flHirano et al. 


19961) . Interferometric observations detected 


a dense gas condens ation toward 


to the outflow axis flOhashi et al 


le pr o 


1996 


osta r with its velocity gradient perpendicular 


19971) . Furthermore Submillimeter Array 


(SMA) observations of this protostellar envelope show that higher-velocity components 
(|1/| 1.4 km s“^) are located closer to the central protostar than the lower-v elocity 


comp onents (|1/| ^ 0.8 km s ^), suggesting a differential rotation of the envelope flYen et al. 


20131) . Their analysis, using the Position-Velocity (PV) diagram with a cut perpendicular 
to the outflow axis to obtain information purely about rotation by avoiding contamination 
from outflowing or infalling gas, showed that TMC-lA has a shallower rotational velocity 
prohle {p = —0.6 ± 0.1 with v oc r^) and higher rotation velocity (1.0 — 4.0 km s“^) than 
their other samples such as L1527 IRS and L1448-mm on a 1000 AU scale. They suggest 
that TMC-lA should have a Keplerian disk {p = —0.5) with M* = 1.1 ± 0.1 Mq and a 
large angular momentum should have been transferred into the disk from the surrounding 
envelope. This result can be understood in terms of the conventional analytical model with 
inside-out collapse. Based on this model, the rotation prohle indicates that TMC-lA is in 
an evoluti o nary stage at t ~ 5 x 10® yr after the onset of collapse. A similar analysis to 


Yen et al. 


( 20131 ) was 


Bure Interferometer (IHarsono et al. 


applied to the visib ilities of TMC-lA data obtained with Plateau de 


20141) . They measured gas kinematics in uv domain 


and suggested that the rotational velocity prohle in the inner part (r < 80 — 100 AU) 
is explained with the power-law index of p = —0.5 (i.e. Keplerian disk) while the outer 
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part is explained with p = —1. They also estimated the size of the disk by removing 
a spherical-infalling envelope component from the continuum emission toward TMC-lA, 
although it was not well constrained (-Rdisk = 80 — 220 AU). These pieces of previous 
work indicate that TMC-lA is in the disk-developing phase, and hence is a good target to 
investigate the formation process of Keplerian disks around protostars. 

Our observations and data reduction are described in SecjH In SecISl we present the 
continuum and molecular-line results. In SecHl we analyze velocity structures of the 
line, in particular rotating motions, while we discuss infalling motions around TMC-lA in 
more detail in SecJSi In SecJHl we present a summary of the results and our interpretation. 
In addition, an investigation into the analysis using a PV diagram and a comparison of 
flattened and spherical envelope models are presented in Appendix and [Bl 


2. ALMA Observations 


We observed our target, TMC-lA during Cycle 0 using the Atacama Large 
Millimeter/Submillimeter Array (ALMA) on November 6, 2012. The configuration is 
“extended” with 23 antennas. Extended conhguration covers projected uv distances from 
16 m to 285 m (11.7—208.7 kA in C^®0 J = 2 — 1). This minimum baseline resolves out 


more than 50% of the flux if a structure is extended more than 777 fjWilner fc Welch 


19M), 


corresponding to ~ 1100 AU at the distance to TMC-lA. The coordinate of map center is 
a(J2000)=4^39“35?010, 5(J2000)=25°4T45"500. ^^CO (J = 2 - 1) and C^^O (J = 2 - 1) 
lines and 1.3-mm continuum emission in Band 6 were observed for 63 minutes (on source). 
In order to derive high velocity resolutions for molecular line observations, we configured 
the correlator in Frequency Division Mode (FDM) for four spectral windows. Each spectral 
window has 3,840 channels covering 234 MHz bandwidth. Emission-free channels are used 
to make the continuum map centered at 225 GHz. The total integrated frequency width 
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of the continuum map is ~ 920 MHz. The observed visibilities were Fourier transformed 
and CLEANed with Common Astronomy Software Applications (CASA). In this process, 
we adopted Briggs weighting with the robust parameter of 0.5 and binned 2 channels thus 
the frequency (velocity) resolution in this paper is 122 kHz corresponding to 0.16 km s“^ 
in ^^CO (J = 2 — 1) and 0.17 km s“^ in C^®0 (J = 2 — 1). We also set CLEAN boxes to 
enclose only positive emission in dirty maps. The synthesized beam sizes of the CLEANed 
maps are 1702 x 0790 for the ^^CO line, 1706 x 0790 for the C^®0 line, and 1701 x 078 7 
for the continuum. J0522-364, J0510+180, and Callisto were observed as Passband, Gain, 
and Flux calibrator. In regard to the flux calibrator Callisto, models of solar system 
objects were updated from Butler-Horizons-2010 to Butler-JPL-2012 for CASA 4.0 after the 
delivery of our data and the fractional difference (2012/2010 —1) was reported. Following 
the value (Callisto: —0.15, 241 GHz), we multiplied the delivered intensities by 0.85. The 
rms noise levels are measured in channels where emission is detected to take the noise due 
to spatial hltering artifacts as well as the thermal noise into account. The parameters of 
our observations mentioned above (and others) are summarized in Table [H 


3. Results 

3.1. Continuum and 

Fig. [T] shows the 225 GHz continuum emission in TMG-lA observed with ALMA. 
Strong continuum emission with a weak extension to the west is detected. Though the size of 
the continuum emission is almost similar to the beam size (1701 x 0787), its deconvolved size 
can be estimated from a 2D Gaussian htting to be 0750 ± 0703 x 0737 ± 0705, P.A.= 73° ±16°. 
The peak position is also measured from the Gaussian htting to be a(2000)=4'^39“35.20®, 
(5(2000)=±25°41'44735, which is consistent with previous measurements (Yen et ah 2013). 
The measured continuum peak position is used as the central protostar position of TMG-lA 


in this paper. The peak intensity and the total flnx density of the continunm emission 
derived from a Gaussian htting are 148.5±2.2 mJybeam”^ and 182.1 ±2.6 mJy respectively. 
The dust mass calculated with the total flux density is Mdust = 
assuming a standard opacity coefficient fi:(l THz) = 10.0 cm ^ g~^ 


opacity index jS = 1.46, and a dust temperature Tc = 28 K fj Chandler et al 


(4.2 ± 0.9) X 10-^ 

(Beckwith et al. 

1990 

), an 


1998|), which is 


a temperature estimated from a htting to the SED of TMC-lA at A 20 fim. Thus, the 
gas mass is estimated to be 4.2 x 10“^ Mq if assuming the gas/dust mass ratio is 100. 

Fig. [2] shows the integrated intensity (moment 0) and the intensity-weighted mean 
velocity (moment 1) maps of the (J = 2 — 1) emission toward TMC-lA, presented 

respectively in contour and in color. The emission is detected at more than 3a level 

at the velocity ranging from —9 km s“^ to 12 km s“^ in Local Standard of Rest frame 
(LSR), except for the velocity from 5 km s“^ to 7 km s“^, where the emission is 

most likely resolved out. This resolved-out velocity range is around the systemic velocity 
hsys = 6.4 km s“^. The emission shows a bipolar structure and a shell-like shape. The 
overall shape and velocity gradient seen in the maps indicate that the emission clearly 
traces a molecular outflow going in the north-south direction. The position angle of the 
outflow axis is measured to be —17° by eye based on the direction of the velocity gradient 
and the shape of the emission, as indic ated with the dashed line in Fig. O which 


is consistent with previous papers (e.g. —10°; 


Chandler et al 


199611 and perpendicular to 


the elongation direction of the continuum emission. The blue lobe is detected more clearly 
than the red lobe. As suggested by the previous ^^CO (J = 1 — 0) observation (Tamura 
et al. 1996), an accelerated motion along the outflow axis can be seen, especially in the 
blue lobe on ~ 1000 AU scale as shown in Fig. [3] where the Position-Velocity diagram of 
the ^^CO emission along the outflow axis is presented. Essentially, the velocity increases 
out along the outflow axis (“Hubble law”) from the central protostar. On the other hand, 
across the outflow axis, the higher-velocity component is located closer to the outflow axis. 
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This feature across the outflow axis has been observed in other outflows as well, and can 


be understood by an outflow driven by a parabolic wide-angle wind (e.g. 


Leeetal 


2000 ). 


In addition to the component tracing the outflow, there seems to be another 

component near the central star, which shows an elongated structure almost perpendicular 
to the outflow axis. This additional component also has a velocity gradient along its 
extension. The geometrical and kinematical structures of this additional component 

are similar to those of (J = 2 — 1), as will be shown later. 


The compactness of the continuum emission with its position angle perpendicular to 
the outflow axis suggests that the continuum emission arises from a compact disk and the 
major axis of the expected disk should be at P.A.= 73° direction (white dashed line in 
Fig. \B). Although the uncertainty of the elongation direction of the continuum emission is 
relatively large (±16°), we hereafter adopt this position angle as the disk major axis around 
TMC-IA. The ratio of minor axis to major axis corresponds to the inclination angle of 
i = 42° and 55° if assuming a geometrically thin disk and H/R = 0.2, respectively. 


3.2. C^®OJ = 2-l 


C^°0 (J = 2 — 1) is detected at more than a 3a level at LSR velocities ranging 
from 2.7 km s“^ to 10.4 km s“^. In Fig. 01 the moment 0 (MOMO) map integrated over 
this velocity range is shown in white contours, overlaid on the moment 1 (MOMl) map 
shown in color. The overall distribution of the C^®0 integrated intensity exhibits an 
elongated structure almost perpendicular to the molecular outflow, with a peak located at 
the protostar position. The deconvolved size of MOMO map is 3.3" ± 0.1" x 2.2" ± 0.1" 
with P.A.= 67° ± 2°. This morphology showing elongation perpendicular to the outflow 
axis (P.A.=—17°) indicates that the emission mainly traces a flattened envelope 


around TMC-lA, as was also suggested by 


Yen et al. 


(1201311 . In addition, there are weak 
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extensions to the north, the northwest, and the east. The overall velocity gradient is seen 
from northeast to sonthwest, which is almost perpendicnlar to the ontflow axis. These 
resnlts are qnite consistent with previons observations in (J = 2 — 1) nsing SMA (Yen 

et al. 2013) althongh their map did not show the weak extensions detected in onr ALMA 
observations. 

In order for ns to see more detailed velocity strnctnres, channel maps shown in 
Fig. [5] are inspected. At high velocities, IAsr < 4.9 km s“^ and > 7.9 km s“^ (which 
are blneshifted and redshifted respectively by more than 1.5 km s“^ with respect to the 
systemic velocity of 6.4 km s“^), there are compact emission with strong peaks located 
nearby the protostar position. The sizes of these emission are smaller than ~ 3" or 420 
AU. The peaks of these emission are located on the east side of the protostar at blneshifted 
velocities, while those for redshifted velocities are located on the west side, making a 
velocity gradient from east to west. This direction is ronghly perpendicnlar to the ontflow 
axis (P.A.= —17°). On the other hand, at lower blneshifted and redshifted velocities 
(5.0 km s“^ < Ylsr < 6.0 km s“^), extended strnctnres become more dominant. At lower 
blneshifted velocities, strnctnres elongated from northwest to sontheast with additional 
extensions to the northeastern side and the sonthern side appear, while at lower redshifted 
velocities, strnctnres elongating to sonthwest become more dominant. At velocities close 
to the systemic velocity, weak emission appears in even more extended strnctnres to east, 
sonth, west, and north, which seem to form an X-shaped strnctnre as a whole. A similar 
X-shaped strnctnre can be also seen in MOMO and MOMl maps althongh they are less 
obvions becanse of a dilntion effect after the channel integration. 

Since the emission shows both compact emission at higher velocities and 

extended emission at lower velocities (whose origin and natnre may be different from each 
other), integrated channel maps of blneshifted and redshifted emission, as presented in 
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Fig. ini were made for further investigation; Fig. |6^ shows channel maps integrated for the 
high-velocity (|A1/| > 2.0 km s“^) compact emission, and Fig. [Hb shows those integrated 
for the low-velocity (|A1/| < 1.0 km s“^) extended emission where AV indicates the relative 
velocity to the systemic velocity. The maps of the high-velocity component show a very 
compact structure as a natural consequence of the integration of compact emission seen in 
the original channel maps shown in Fig. 0 The dashed line in Fig. passing the two peaks 
of blue- and redshifted components gives us a direction of the velocity gradient at the high 
velocities (P.A.= 72°), which is almost perpendicular to the outflow direction. In contrast 
to the high-velocity component, the maps of the low-velocity component are extended to 
~ 800 AU scale at the Aa level. This hgure also shows a clear velocity gradient in the 
direction of northeast-southwest as shown with a dashed line (P.A.= 49°) passing the two 
peaks of blue- and redshifted components in Fig. Eb. This direction of the velocity gradient 
of the low-velocity component is not perpendicular to the outflow axis and apparently 
different from that of the velocity gradient seen in the high-velocity component. These 
differences seen between the high-velocity and low-velocity components will be discussed in 
more detail later. 

The ratio of the brightness temperature between the and emission tells 

us that is most likely optically thin although it is difficult for us to estimate the 

optical depth of the emission at between 5 km s“^ and 7 km s“^, where is 

completely resolved-out. We note that when the ratio between the and emission 
is calculated, the spatial hltering has to be taken into account. In the case discussed here, 
however, the emission is signihcantly more extended than the emission, and 

as a result, the emission should be resolved out more than the emission. This 

suggests that the ratio of the emission to the emission should be overestimated, 
and as a result, the optical depth of the emission should be overestimated. When 

is optically thin and under the condition of Local Thermodynamic Equilibrium (LTE), 
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total gas mass is estimated from total flux of C^®0. The peak intensity and the total flux of 
the emission are 842 ± 82 mJy beam“^ km s“^ and 7.20 ± 0.70 Jy km s“^, respectively. 
The gas mass calculated with the total flux is Mg^s = 4.4 x 10“^ Mq assuming the excitation 


temperature T = 28 K (Chandler et a. 

= 3.0 X 10“^ ( Frerking et ah 


■ 1998 ) and the abundance of relative to H 2 of 
19821 ). The gas mass derived from the flux is 


one order of magnitude smaller than that from dust conti nuum flux density^ 4.2 x 10 ^ Mq 


assuming gas/dust=100 (see a result for Class 11 disks by 


Williams fc Best 


2014j) . though 


both mass estimations include a lot of uncertain factors such as optical depth, temperature, 
non-LTE effects, and gas/dust mass ratio. Additionally the molecular abundance of C^®0 
might be decreased because of freeze-out of CO onto dust grains. 


4. Analysis 

4.1. Velocity Gradient 

As shown in the previous section, the C^®0 emission arising mainly from the flattened 
envelope of the TMC-lA protostar shows a very clear velocity gradient. In this section, the 
nature of the velocity gradient is investigated in more detail. 

One remarkable characteristic of the velocity gradient is that its direction changes 
between high and low velocities as shown in Fig. [6l In order for us to investigate such a 
tendency more systematically (channel by channel), the mean position of the C^®0 emission 
at each channel is measured, and plotted in Fig. [71 The mean position is measured as 
*o('i^) = / xl{x,v)(fx/ f I{x,v)(fx, where the sum is calculated with pixels having an 
intensity more than 6cr [a is the rms noise level of C^^O). For channel at |AV| < 0.5 km s“^, 
the mean positions are not measured because very complex structures appear at those 
channels, such as the “X shape”. Points corresponding to higher velocities are plotted with 
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smaller point sizes in Fig. [7] and green points correspond to |A1/| > 2.4 km s“^. This 
figure demonstrates the difference between high and low velocities seen in Fig. [Hlmore 
systematically, i.e., the mean positions for higher velocities are mostly along the major axis 
of the disk (dashed line, E-W) and near the center. Those for lower velocities depart from 
the major axis and are located far from the center. This displacement changes the direction 
of the velocity gradient between higher and lower velocities. Fig. [7] suggests that the motion 
in the inner (high-velocity) part is dominated by rotation of the disk because the velocity 
gradient is almost along the major axis of the disk. In contrast to the motion in the inner 
part, the motion in the outer (low-velocity) part is considered to include additional motions 
to rotation, which give a tilt in the velocity gradient with respect to the disk major axis. 

In order for us to examine velocity gradients along the disk major axis in more detail, 
the Position-Velocity (PV) diagram cutting along the disk major axis shown in Fig. is 
inspected here. This figure shows that the velocity gets higher as the position gets closer 
to the central star. As we discussed above, the PV diagram along the disk major axis is 
naturally considered to represent a velocity gradient due to rotation of the disk and/or the 
disk-like flattened envelope. There is a feature of the so-called “spin up” rotation where 
emission peaks get closer to the central position (horizontal dashed line) as the velocity 
increases, while emission peaks depart from the center at low velocities. Keplerian rotation 
(Vot cc is one form of spin up rotation, and rotation that conserves its specific 

angular momentum (Vot oc R~^) also shows a similar spin-up feature. However, it is 
generally not easy to distinguish such rotation patterns by comparing them with the PV 
diagram along the disk major axis. Further analysis to determine the powers of rotational 
velocity will be presented in the following section. 
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4.2. Rotation Profile 


In order for us to characterize the nature of the rotation, detailed analysis of the PV 
diagram is presented in this section. The nature of the rotation can be characterized by its 
rotation prohle obtained from the PV diagram by measuring a representative position at 
each velocity channel in the diagram. The representative position at each velocity channel is 
dehned as ID mean position weighted by the intensity i.e. Xmiv) = / xl{x, v)dx/ f I(x, v)dx 
where the sum is done with pixels having an intensity more than 6a (a is the rms noise 
level of C^®0). The measured mean positions are superposed on Fig. [8^ with error bars. 
The error bar along the direction of the position is derived by the htting while one along 
the direction of the velocity is the velocity resolution of the C^®0 observations. The mean 
position is not measured at low velocities (|A1/| < 1.2 km s“^) where the emission shows 
extended structures which is unsuitable for this analysis. We also note that a similar 


analysis was performed by 


Yen et ah 


(120131). 


Fig. |9] shows the mean positions of the PV diagram along the major axis on a 
log R — log V diagram. The abscissa shows the mean positions measured as offset positions 
from the central protostar, which corresponds to the zero position on the PV diagram, 
assuming that the distance of TMC-lA is 140 pc. The ordinate is the relative-projected 
velocity assuming that the systemic velocity is 6.4 km s“^. Because there is no systematic 
separation between the blue- and red-shifted components, we consider that the central 
position and the systemic velocity are plausible in order to measure the relative position and 
velocity. Fig. [9] shows clear negative correlation between radius and velocity, demonstrating 
that the rotation velocity can be expressed as a power-law function of the radius (spin-up 
rotation). Importantly, there seems to be a break point at i? ~ 70 — 80 AU, where the 
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power changes. We performed a fitting to Fig. |9]with a double power model as follows: 

[ (R > Rt) 

The powers, Pin and Pout, and the break point {Rh, 14) are free parameters. The best fit 
parameter set is {Rb, 14, Pin, Pout) = (67 AU, 2.4 km s“^, 0.54 ± 0.14, 0.85 ± 0.04) where 
the errors for Rb and 14 are ~ 1%, giving reduced = 1-6- Because this reduced is 
close to the unity, we consider that the best £t is reasonable. The error of each parameter 
is derived by dispersing the data points randomly within the error bars of logi? and log 17. 
We adopted a Gaussian distribution for the error of log R and a rectangular distribution 
for the error of log 17 because the error of log 17 is based on the velocity resolution of the 
observations. Another htting to the mean positions is performed with a single power law. 
The best-fit parameter set for this single power law is {Rb, p) = (64 AU, 0.69 ±0.02) giving 
reduced = 2.0. In the single-power case, Vb is hxed at 2.4 km s“^ and {Rb, 14) does not 
give a break. Although more parameters usually make a htting better or x^ smaller, in our 
case, the best £t with a double power law provides a smaller “reduced” than that with a 
single power law even taking account the difference in the number of parameters between 
two httings (four for the double power law and two for the single power law). This suggests 
that the rotation profile of TMC-lA can be htted with the double power law better than 
the single power law. The best-£t result suggests that the inner part and the outer part 
follow different rotation laws. The power of the inner part is roughly equal to 1/2 for the 
case of Keplerian rotation, suggesting that the inner/higher-velocity part is dominated by 
Keplerian rotation. By contrast, the outer/lower-velocity part, showing a slope steeper 
than Keplerian, cannot be supported by rotation against the gravity yielded by the central 
protostar, as we will discuss later. It is interesting to note that Pout is equal to —0.85, 
which is signihcantly different from —1 for the case where the specific angular momentum 
is conserved. This result could suggest that the specific angular momentum is proportional 
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to 

By using the break point {Rb, H) = (67 AU, 2.4 km s“^), the dynamical 
mass of the central protostar and the specihc angular momentum at Rb, which are 
not inclination corrected, can be calculated at = M* sin^ i = 0.43 Mq and 
jp = j sini = 7.7 x 10““^ km s“^ pc, respectively. We should note, however, that the analysis 
using certain representative points on a PV diagram presented here may systematically 
underestimate Rb when the spatial resolution is not high enough, as discussed in Appendix 
lAl Taking this point into consideration, the real break radius can be estimated to be ~ 90 
AU from the apparent break radius of 67 AU. Rb ~ 90 AU leads ~ 0.6 Mq and 
jp ~ 1 X 10“^ km s“^ pc. 


4.3. Rotating Disk Models 


In the previous section, the detailed analysis of the PV diagram cutting along the 
disk major axis was presented to characterize the nature of rotation, which suggests the 
existence of a Keplerian disk in the innermost part of the envelope. In this section, physical 
conditions of the Keplerian disk are investigated based on model httings to the 
channel maps. To quantify physical parameters of the disk around TMC-lA, we performed 
a htting to the J = 2 — 1 channel maps based on a standard rotating disk model 


(e.g. 


Dutrev et al. 


I994J) . When a standard disk model is compared with these observed 


channel maps, it is important to note that the observed C^®0 emission arises not 
only from the disk, but also from the envelope having rotation and additional motions as 
discussed in Sec J4.ll Because of this, it is required that only the velocity channels where 
rotation is dominant should be used for comparison with a disk model. As shown in Fig. [HI 
there is a velocity range showing the power-law index close to the Keplerian law, suggesting 
that the emission in this range with |AU| > 2.4 km s“^ is thought to arise from the 
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possible Keplerian disk. We, therefore, intend to fit the channel maps within this velocity 
range, with rotating disk models. 


We use a disk model with a code described in lOhashi et al.l (120141) . The model disk 
can be described with 10 parameters summarized in Table [2l The radial dependence of the 
disk temperature T{R) and the disk surface density S(i?) are described as 


T{R) = T] 


100 


R 


100 AU 


-9 




- Kn 

where M 200 is the mass within 200 AU (i?out will be hxed at 200 AU, see below for more 
detail). The scale height of the disk H{R) is given under the hydrostatic equilibrium 
condition as follows: 


( 2 ) 

( 3 ) 


H{R) = 


2kBT{R)R^ 


( 4 ) 


Gm^M^p/ sin^ i 

where ks, G, uiq, and i are the Boltzmann constant, gravitational constant, mean molecular 
weight (=2.37 mu), and inclination angle of the disk, respectively. In such a disk, its density 
distribution p{R,z) is given as S(i?)/(i/^i7(i7)) exp(—2;^/i7(i?)^), while the rotation prohle 
of the disk is provided by 


VtiR) = 


' GM^p/ sin^ i 


R 


( 5 ) 


Rccnt V -^cent 

where i?cent is set to satisfy V^(i?cent) = GM^,p/ sin^ i/Rcent and gives a coefficient of the 


power law of V^{R). It is assumed that the turbulent linewidth is zero; that is, the prohle 
function is dehned by a thermal linewidth. 


Among the 10 parameters, the hrst three parameters {i, R,^, i?out) are hxed, as shown 
in Table Hi Note that i?out=200 AU is assumed to cover a sufficiently large area for the 
htting because our htting to only the high-velocity inner part cannot constrain Ront well. 
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In addition, i = 65° is also snitable to explain motions in the envelope as discussed later 
in Sec J5.2[ With the other seven free parameters {M^p, M 2 oo,p,TiQQ,q, Rcent,v), our model 
solves radiative transfers in 3D + velocity space to produce model channel maps (data 
cube). When solving the line transfers, we assume Local Thermodynamic Equilibrium 
(LTE) and a typical abundance of relative to H 2 , X(C^®0) = 3.0 x 10“^. After solving 
the line transfers, model data cubes are convolved with the Gaussian beam, which has the 
same major axis, minor axis, and orientation as the synthesized beam of our observations. 
When the model is compared with the observations, the central position and the orientation 
(position angle, P.A.) of the disk are assumed to be the same as those of the continuum 
emission and the systemic velocity is also assumed to be 14ys = 6.4 km s“^. 

Reduced is calculated to evaluate the validity of each model in the velocity range, 
|AR| > 2.4 km s We use only pixels in channel maps (data cube) where observed 
emission is stronger than the 3a level because, if undetected pixels are included, where noise 
is distributed around la (reduced ~ 1), then the value would be underestimated. We 
use Markov Chain Monte Carlo (MCMC) method to hnd the minimum y^ efficiently. 

Fig. ITOh shows the comparison of the best-ht disk model with the observations and 
Fig. ITOb shows the residual obtained by subtracting the best-ht disk model from the 
observations. Reduced y^ = 0.72 and Fig. flOb indicate that our best-ht disk model 
reproduces, overall, the observations in the high-velocity region. The parameters of our 
best-ht disk model are summarized in Table [2l First of all, v = 0.53 suggests that the 
high-velocity component of C^®0 emission from TMC-lA can be explained as a Keplerian 
disk better than other rotation laws, such as oc R~^. Note that this power-law index is 
very consistent with pin derived from the power-law htting to the log R — log V diagram in 
Sec J4.21 The central stellar mass (not inclination corrected) = M*sin^f = 0.56 Mq is 
also consistent with that derived in Sec J4.2[ With the best-ht M*p, v, and i?cent = 166 AU, 
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the rotational velocity of our best-fit model is expressed as 1.5 km s“^(i?/200 The 

mass within i?out = 200 AU, 3.6 x 10“^ Mq is similar to the total gas mass derived from 
the flux of C^®0, Mgas = 4.4 x 10“^ Mq. With the best-£t M 200 and p = 1.46, the surface 
density prohle of our best-£t disk model is expressed as 0.069 g cm“^(i?/200 AU)“^'^®. The 
temperature seems roughly uniform [q ~ 0) in our best-fit disk model. The temperature at 
the midplane can show such a radial prohle when the inner region has a larger optical depth 
enough to cause a more eff ective cooling and a less effective heating by hot surface regions 


fChiang fc Goldreic 




19971). Another possibility is that the line traces a temperature 


at which C^®0 molecules are evaporated off of dust grains. 


The uncertainty of each parameter is dehned by the range of the parameter where 
the reduced is below the minimum plus one (= 1.72) when varying the parameter and 
hxing other parameters at those of the best-ht disk model. M* and v, which are related 
to kinematics, have smaller uncertainty as compared to others. Another parameter related 
to kinematics, Rcent, has a large uncertainty because the rotational velocity is found to be 
close to the Keplerian rotation law. 

It is important to investigate the effect of the inclination angle, which is hxed at 
i = 65° in our htting mentioned earlier. Fig. [TT] shows a distribution of the reduced 
with different inclination angles of the disk. This plot is derived by changing the inclination 
angle i from 5° to 85° by 5° intervals, with the other parameters hxed at those of the 
best-ht disk model. Fig. [TT] shows that larger inclination angles between 45° and 75° give 
small reduced y^, suggesting that it is difficult to give a constraint on the inclination angle 
in this model htting. In Sec. 15.21 we examine the disk inclination angle further. 
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5. Discussion 

5.1. Velocity gradient along the disk minor axis 

In the previous sections, detailed analysis of the motions for the inner part of the 
envelope where rotation is dominant are presented, showing that there is a Keplerian 
disk at the innermost envelope. As shown in Sec. 14.11 the kinematics of the outer part of 
the envelope surrounding the Keplerian disk is described with not only rotation but also 
additional motions. In this section, the nature of the additional motions in the outer part 
of the envelope is discussed in detail. 

There are two possibilities for the additional motions seen in the outer part. One is 
radial motion in the disk plane and the other one is motion perpendicular to the disk plane. 
Either case is expected to form velocity gradients along the disk minor axis. In the latter 
case, outflowing motion going perpendicularly to the disk plane is the most natural one. 
On the other hand, in the former case, infalling motion or expanding motion in the disk 
plane can be considered, but these two can be distinguished from each other because the 
directions of their velocity gradients are opposite to each other. The geometry of the outflow 
having its blueshifted lobe located on the northern side of the central star tells us that the 
northern part of the disk plane is on the far side from the observers, while the southern 
part is in the near side. Since the northern part of the emission is expected to have 

blueshifted emission as mentioned in Sec. 13.21 the radial motion in the disk plane should be 
infall rather than expansion. Based on these considerations, the additional motion would 
be either an infall motion in the disk plane, an outflow motion going perpendicularly to the 
disk plane, or both. 

The PV diagram along the disk minor axis shown in Fig. [Hb shows a velocity gradient 
with blueshifted and redshifted components located on the northern side and the southern 
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side, respectively at low velocities (4.0 km s“^ ^ 14 8.8 km s“^). On the other hand, 

high-velocity {V ^ 4.0 km s“^, 8.8 km s“^ ^ 14) components show no signihcant velocity 
gradient along the disk minor axis, which is consistent with what was discussed with Fig. [3 
The velocity gradient seen at low velocities can be explained as either infall motion in the 
disk plane, outflow motion going perpendicularly to the disk plane, or both, as discussed 
above. It should be noted that a molecular outflow needs velocities sufficiently high enough 
to escape from the gravitational potential due to a central star. Such escape velocities are 
estimated as a function of the distance from the central star and are potted in Fig. |8 ]d. 
When plotting the escape-velocity curves at the low velocities where there is a velocity 
gradient, the central mass (not inclination corrected) = M*sin^i = 0.56 Mq was 
assumed based on the estimation of the dynamical central stellar mass (see Sec. 14.31). T he 


inclination angle was also assumed to be 40° < i < 70° based on 


Chandler et ah 


(119961). 


Compared with the curve with i = 70°, a part of the emission shows higher velocities 
than the escape velocities, suggesting an outflow motion. In fact, a part of the emission 
shown in channel maps (Fig. [5]) exhibits extensions or X-shaped structures, which spatially 
correspond to the molecular outflow traced by ^^CO. The signihcant part of the emission 
including the two peaks in Fig. |8 ]d, however, show clearly slower velocities than the escape 
velocities. This suggests that the major part of the velocity gradient along the minor axis 
cannot be explained by outhow motions. We therefore consider that the velocity gradient 
along the disk minor axis in C^®0 emission is mainly due to infall motion in the disk plane. 
Note that additional signature for the infall motion can be also seen in the major-axis 
PV diagram at lower velocities (5.5 km s“^ ^ C ^ 7.0 km s“^): although the most of the 
blueshifted emission is located on the eastern side, a part of the blueshifted emission at 
lower velocities is located on the western side, and similarly a part of the redshifted emission 
at lower velocities is located on the side opposite to the western side where the most of 
the redshifted emission is seen. These additional components in the diagram cannot be 
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explained by rotation, and can be explained by infall (e.g. 


Ohashi et al. 


1997 1. In addition, 


we also note that L1551 IRS5, a typical Class I pr otostar, also shows a n infalling envelope 


with a velocity gradient along its disk minor aixs flMomose et al 
snbsections, details of infall motions seen in are discussed. 


19981) . In the following 


5.2. Nature of the infall motion 

It is suggested that the inner high-velocity part of the emission arises from a 

Keplerian disk and the parameters of the disk are derived by model httings in Sec. 14.31 
In addition to the rotation-dominant part, our results suggest that there is an infalling 
envelope surrounding the rotating disk. In order to investigate infall motions in the 
envelope, we make another model by adding an infalling envelope to the standard rotating 
disk model used in Sec. 14.31 Although the disk is basically described based in the best-fit 
disk model obtained in Sec. 14.31 it is necessary to define the transition radius between the 
disk and the infalling rotating envelope because the outer radius in the best-£t disk model 
is simply assumed to be 200 AU. A radius i?kep is set as a boundary between the infalling 
rotating envelope and the purely rotating disk. Within Rkep, all the parameters are the 
same as the best-fit disk model obtained in Sec. 14.31 except for the inclination angle (i) 
which was not fully determined in Sec. 14.31 The velocity held of outer part {R > i?kep) is 
described as below: 

vm = 

V(r) = 

Note that we distinguish the axial radius R = and the spherical radius 

r = by using big and small letters, respectively. The infall-velocity vector 

points to the center and the infall velocity Vr{r) depends on z while the rotation velocity 


h^(77kep) 


R 


R 


kep 


a 


2GM^p/ sin^ i 


( 6 ) 

(7) 
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V^{R) does not depend on it. V"^(i?kep) is the rotation velocity at i?kep and is assnmed 
to connect continnously to that of the inner disk. The infall velocity is set as a product 
of a constant coefficient a and the free fall velocity. In terms of the geometry of the 
infalling envelope, we basically consider a flattened structure because of the shape of the 
integrated-intensity map (Fig. |4j). This flattened geometry of the infalling envelope 
is modeled by extrapolating the density and temperature structures from the inner disk 
{R < i?kep) to the outer radius of 900 AU, which is larger than the size of the 
integrated-intensity map. Hereunder, the models with a flattened envelope are called FE 
model. 

It is important to note that in Sec. 14.31 the rotating disk models are htted to the 
channel maps at higher velocities, while such httings are not performed here. This is 
because models including the infalling envelope must be compared with the channel maps 
not only at higher velocities, but also at lower velocities with complicated asymmetric 
structures. Such structures cannot be fitted by our simple axisymmetric models. For this 
reason, in this section, we examine whether our models can reproduce the observed PV 
diagrams cutting along the major and minor axes (Fig. [8]), by comparing the following two 
parameters between them: (1) the mean position in the major-axis PV diagram measured 
in Section 02] (Fig. [Hti' and Fig. [9]) and (2) the position of two peaks in the minor-axis PV 
diagram (Fig. |8 ]d). The first depends on the rotation velocity while the second depends on 
the infall velocity within the disk plane. 

It is also noted that the surface density of the envelope at radii more than 200 AU is 
reduced by a factor of two to reproduce relatively diffuse emission seen in PV diagrams 
at positions larger than ~ 2". Because extended structures more than ~ 8" are resolved 
out in our ALMA observations by 50%, we use the CASA task simobserve to include the 
spatial hltering. The synthesized beam derived with simobserve (1'.'03 x Cf.'81, —7.2°) is not 
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exactly the same as that of the observations even though the same antenna positions and 
hour angles are conhgured. However, the difference is small and there is no impact on our 
analysis. 

Four parameters of the model mentioned above, Uout, .Rkep, «, and the inclination 
angle of the disk i are changed to reproduce the observed PV diagrams. Fig. [T2] shows 
the comparison of different models (FE model A-F) and the observations. The FE model 
A is the case with (uout, .Rkep, a, *) = (0.85, 100 AU, 0.3, 65°). This FE model A can 
reproduce both PV diagrams, whereas the other hve models cannot reproduce either of 
them as explained in detail below. Hereafter, we call the FE model A “best-infall model”. 
We emphasize that the infall velocity is signihcantly smaller than the free fall velocity in 
the best-infall model. More details of the slow infall velocity is discussed later. 

FE model B and C, which have a shallower Uout and a smaller i?kep 5 respectively, as 
compared with FE model A, show the clear difference in the log R — log V diagram. The 
mean position derived from FE model B are located above the observations in i? 70 AU. 
On the other hand, the mean positions of FE model C clearly show smaller radii than the 
observations in i? 50 AU. 

FE model B and C exhibit similar PV diagrams along the minor axis to that of FE 
model A. However, FE model D and E, which have a larger infall velocity and a smaller 
inclination angle, respectively, as compared with FE model A, show clear difference in 
the minor-axis PV diagram. The two peaks of FE model D shows a clear offset from the 
observations along the “velocity” axis in the PV diagram. In addition, the log R — log V 
diagram of FE model D is signihcantly different from the observations, in the sense that 
velocities at radii less than ~ 80 AU become larger than the observations, which can be also 
seen in the major-axis PV diagram where the two peaks in the model have larger velocities 
than the observations. On the other hand, the two peaks in the minor-axis PV diagram of 
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FE model E shows a clear offset along the “position” axis in the PV diagram. 

The fact that FE model B, where the outer and inner power-law indices are the same as 
0.53, cannot reproduce the observations indicates that the power-law index of the rotational 
velocity in the outer region is steeper than that in the inner region. The discrepancy 
between the observations and the FE model C with i?kep = 70 AU, at which the observed 
rotation profile appears to break, indicates that i?kep should be larger than the apparent 
breaking radius, as also discussed in Appendix (see also Sec. 14.2p . The discrepancy 
between the observations and FE model D, where the infall velocity is free fall velocity, 
indicates that the infall velocity in the envelope should be smaller than the free fall velocity. 
The fact that the FE model E with i = 50° cannot reproduce the observations indicates 
that the inclination angle of the disk is larger than i = 50°. Though the inclination angle is 
not well constrained with the disk model, as shown in Sec. 14.31 it is constrained better in 
FE models because the peak positions in the minor-axis PV diagrams are sensitive to the 
inclination angle. 

In addition, the power law dependence of the infall velocity is investigated. Although, 
in FE model A, a conventional infall velocity law with the power-law index —0.5 is assumed, 
it would be possible that infall motions are decelerated and cease at the boundary radius 
because the inner region shows Keplerian rotation and no infall motion. Thus, we make FE 
model F with the same parameters as FE model A except for the infall velocity, which is 
linearly decelerated K-(r) = a(r — i?kep) (-R > -Rkep) where a is a constant coefficient and set 
to be a = 1.6 x 10“^ km s“^ AU~^. Overall velocity structures seen in both PV diagrams 
made from FE model F seem to be consistent with the observations, even though the peaks 
in the minor-axis PV diagram of this FE model F are located at a smaller velocity than 
the observations and the two peaks are almost merged. This would suggest that the infall 
velocity may decelerate, although such a deceleration can happen only at the innermost 
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region of an envelope. 

From the parameters of the best-infall model i.e. FE model A, (fout, -Rkep, *) = 
(0.85, 100 AU, 0.3, 65°), the infall velocity and mass infall rate M can be estimated. 
First, the inclination-corrected central stellar mass is calculated at M* = 0.68 Mq from 
the inclination i = 65° and the best-£t = M^sin^i = 0.56 Mq (best-ht disk model in 
Sec i4.3p . Thus, from a = 0.3, the infall velocity is estimated at K = 0.74 — 1.04 km s“^ at 
R = 100 — 200 AU. Using the surface density of the best-infall model, the mass infall rate 
is estimated as 


M = 2'KRVrT. = (1.5 - 3.0) x 10“^ Mq yr“^ {R = 100 - 200 AU). 


( 8 ) 


Ohashi et ah 


1997h 


This value is comparable with a typical mass infall rate of protostars (e.g. 
and the mass flow rate Maow = 1.5 x 10“® Mq yr“^ of TMC-lA estimated by an observation 
of the molecular outflow (Chandler et ah 1996). The inclination-corrected central stellar 
mass M* = 0.68 Mq and the boundary radius i?kep = 100 AU also give us a specihc angular 
momentum at the outer radius of the Keplerian disk, j = 1.2 x 10“^ km s“^ pc, which is 


consistent with 


a factor of two fj Ohashi et ah 


he specific angular momentum at 580 AU, j = 2.5 x 10 km s ^ pc within 


19971 ). From i?kep = 100 AU, the disk mass can be calculated 


to be Mdisk = / Ti27rRdR = 2.5 x 10 ^ Mq, as well. 


One might wonder whether a spherical geometry of the infalling envelope could explain 
the observations better, including the slower infalling velocities, even though the envelope 
has a flattened structure as shown in Sec. 13.21 We therefore consider models with a 
spherical envelope having bipolar cavities in addition to FE models discussed above as 
shown in Appendix [Bl It is found that models with a spherical envelope cannot explain 
the observations better than EE models. More importantly even if we adopt a model 
with a spherical envelope, the infall velocity should be signihcantly slower than the free 
fall velocity. In addition, it might be possible that power-law index q for the temperature 












27 


distribution in an envelope is not the same as that in a disk. A smaller value of q, i.e., 
—0.5 is examined in the model and found that it does not signihcantly change the results 
discussed above. 


5.3. Magnetic Pressure and Tension 


In previous sections, we suggest that the infall velocity in the envelope around TMC-lA 
is ~ 0.3 times as large as the free fall velocity, Vg, yielded by the stellar mass, which is 
estimated from the Keplerian rotation. The slow infall velocity suggests a possibility that 
infalling material are supported by some kind of force against the pull of gravity. One 
possible mechanism is the magnetic effect. Dense cores where protostars form tend to have 
obl ate shapes perpend icular to their associated magnetic helds. (e.g. NGC 1333 IRAS 


4A, 


Girart et ah 


20061) . Thus, it is assumed that envelopes around protostars should also 


be penetrated by magnetic helds, and the magnetic force (pressure and tension) possibly 
slows infall motions. In this section, we discuss whether the magnetic held expected to be 
associated with the envelope around TMG-lA can actually make the infall velocity as small 
as 0.3 Vff. Note that all the physical quantities in the following discussion are derived at a 
radius of 200 AU, which is rather arbitrary. The same estimations are also done at radii of 
100 AU and 300 AU, providing results similar to that at 200 AU. 

The magnetic hux density required to make the infall velocity as small as 0.3Ur is 
estimated from the equation of motion including magnetic helds described below under a 
symmetric condition of d/d6 = 0 and d/dz = 0 on the midplane: 


p(R}K(R)^ = 


GM.p(R) 

R2 


1 dB^ 


+ 


B{Rf 


2^0 dR /^O-^curv (-R) 


( 9 ) 


where p(R), Vr{R), and i?curv(-R) indicate radial dependence of the mass density and the 
infall velocity in the midplane of the envelope, whereas i?curv(-R) is the radial dependence of 
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the curvature of magnetic field lines, /tq indicates the magnetic permeability in a vacuum. 
The second and third terms of the right hand side indicate the magnetic pressure and the 
magnetic tension, respectively. A steady state is assumed here and the gas and radiation 
pressnres are ignored because they are too weak to support material against gravity, i.e., 
the thermal energy is smaller than the gravitational potential by orders of magnitude at 
~ 200 AU). 


Although we do not know the conhguration of magnetic helds around TMC-lA, it is 
reasonable for our order estimation to assume that Rcuw is roughly equal to the scale height 
H, which is a typical spatial scale in the disk-like envelope. We also assume B oc 
which corresponds to the case where gases spherically symmetrically collapse with magnetic 
helds. Then, with p{R) and Vr{R) derived from the best-infall model discussed in the 
previous section, the magnetic hux density required for W = O.Sldr is estimated to be i? ~ 2 
mG. The magnetic hux density on the order of mG is somewhat larger than the typical 
value from Zeeman observations 100 pG). However, the Zee man observations usi ng a 


high-density tracer snch as GN give us higher values, ^ 1 mG f Grutcher et a 


Girart et ah 


201011 . In 


(1200611 toward 


addition to single-dish observations, observations using SMA by 
NGG 1333 IRAS 4A at ~ 1" angular resolution estimated the magnetic held strength in 
the plane of the sky to be 
These observations (see also 


pos ~ 5 mG bv applvine the Ghandrasekhar-Fermi equation. 

Falearone et al. 

OC 

o 

o 

CM 

Stephens et al. 

2013 

) suggest that the 


magnetic held strength around TMG-lA can be possibly on the order of mG. 


5.4. A Possible scenario of the Keplerian Disk Formation 

Over the last decade the number of protostellar sources associated with Keplerian disks 
have been increased with (sub)milhmeter interferometric observations (see Table 3), and it 
is getting more feasible for us to discuss evolutionary processes of Keplerian disks around 
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protostars based on these observations (e.g. IChou et al.ll2014l) . We should, however, note 
that disk radii measured from such observations are subject to serious ambiguity, because 
the measured radii are likely affected by the effects of the sensitivity limits and the missing 
fluxes of the interferometric observations. In particular, some of these observations did not 
identify the transition from infalling motions to Kepler motions (see Column 8 in Table 3), 
so that the estimated disk radii measured from such observations should be lower limits 
to the actual disk radii (see more details below). Indeed, 


Chou et al. 


fl2014j) attempted to 


discuss evolutionary processes of Keplerian disks around protostellar sources with those 
samples, but they could not identify any clear correlations and evolutionary sequences of 
the disks. 


by 


On the othe r hand , our series of ALMA observations of protostellar s ources, L1527 IR S 


Ohashi et al. 


( 2014 ). TMC-lA in the present work, and L1489 IRS by 


Yen et al. 




all identihed the transitions from infalling envelopes to disks, and measured the power-law 
indices of the rotational prohles to verify that the disks are indeed in Keplerian rotation. In 
such observations, the disk radii are estim ated better. Furtherm ore, recent hig h-resolution 


interf erometric observat ions of L1551 NE flTakakuwa et al. 


2014 ). and VLA1623A ( Murillo et al. 


2014 ). L1551 IRS 5 f Chou et al. 


20131) have also revealed transitions from infalling 


motions to Keplerian rotations. These observational results imply that transitions from 
infall motions to Keplerian rotations are ubiquitous in the protostellar phase. Statistical 
studies of such systems should provide important observational constraints on the ongoing 
growth process of Keplerian disks. Here, with our ALMA results, along with the other 
observational results showing transitions from infalling envelopes to Keplerian disks, we 
re-examine evolutionary processes of Keplerian disks around protostellar sources embedded 
in infalling envelopes. 


In the following discussion, the bolometric temperature {T^oi) and the ratio of the 
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bolometric to submillimeter l uminosity (= Lhni ! are adopted as parameters tracing 


the evolution of protostars. flGreen et al. 


2 OI 3 I ) Fig. [15] shows correlations between T\ 


bol 


VS the central protostellar masses (= M*), L^oi/Lsnhmra vs M*, and M* vs the radii of the 
Keplerian disks (= i?kep) for the protostellar sources associated with the Keplerian disks 
listed in Table 3. Note that the protostellar sources “with transitions” showing transitions 
from infall to Keplerian and “without transitions” are indicated with different marks. The 
sources with transitions exhibit correlations in all the T^oi vs M*, Lboi/Lsuhram vs M* and 
M* vs i?kep plots, with the correlation coefficients of 0.97, 0.95, and 0.91, respectively. The 
correlation coefficient for M* vs i?kep will be 0.98 if assuming that VL A1623A has = 50 


AU, at which its rotational velocity prohle shows a tnrn over point (IMurillo et al. 


201311 . 


These resnlts may snggest that the central protostellar masses increase as the protostellar 
evolntion proceeds, and that the disk radii also grow as the central protostellar masses 
grow, although statistical signihcance are not high because the number of sample is small. 
On the other hand, inclusion of all the sample points almost diminishes the correlation in 
M* vs i?kep with the correlation coefficient of 0.29. This poor correlation may be because 
i?kep is not measured accurately for the cases without transitions. The correlation between 
Lboi!Lsnhrara and M* is not diminished (the correlation coefficient is 0.77) because measuring 
M* does not depend on whether there is a transition thongh the correlation of T^oi vs M* is 
poor when all the samples are included (the correlation coefficient is 0.43). These results 
demonstrate that proper observational estimates of the central protostellar masses and the 
disk radii are essential to discuss evolutions of the protostars and disks. 

A simple interpretation of the positive correlation between M* and i?kep is given by a 
conventional analytical picture of disk formation with “inside-out collapse” of envelopes. 


flTerebev et al 


1984; 

Shu et al. 

1987; 

Yen et al. 

2Q13) 


201311 . At the onset of the collapse, only the 


central part of the envelope with less angular momentum falls toward the center, and thus 
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the mass of the protostellar source and the radius of the Keplerian disk are smaller. As the 
collapse proceeds outward, more material with larger angular momenta will be carried into 
the central region, raising the protostellar mass and radius of the Keplerian disk . In this 
pictur e, the dependence of i?kep on t (time) and M* can be expressed as shown in 

(bold): 


Yen et ah 


2^4 


-^kep 




inf 




t^ = 


-Mt 


( 10 ) 


GM*(rinf) 16cf 16cf 

where oj and Cg denote the angular velocity and the sound speed of the envelope, 
respectively, rinf is the outermost radius of the infalling region, and M is the mass infalling 
rate. A typical yalue of u has been estimated as cu ~ 1 km s“^ pc“^ = 3.2 x s“^ 


((Goodman et al. 


1993). In Fig. [T3b lines of Eq. llTOl) with twice and half of the typical 


value of uj are drawn, with Cg = 0.19 km s“^ (i.e., = 10 K). Almost all the sources with 

transitions are within the area expected from Eq. [TOl suggesting that the identihed positive 
correlation between M* and i?kep can be reproduced with the picture of the conventional 
analytical disk-formation model with inside-out collapse. 

More realistically, effects of magnetic helds on disk formation a nd growth, which ar e 


ignored in the model described above, must be taken into account ((Machida et al. 


Li et al. 


2914 : 


2014 1. Our ALMA observations of TMC-lA and L1527 IRS have found that the 


infalling velocities are a factor of a few slower than the corresponding free fall velocities 
derived from the central protostellar masses, while in L1489 IRS the infalling velocity is 
consistent with the free fall velocity. As discussed in the previous section, the slow infall 
velocities might be due to magnetic helds, while in the case of L1489 IRS, the magnetic 


held may not be ehective any more because the surrounding enve 


magnetic held is almost dissipat ed in the late evolutionary s 


latest theoretical simulations by 


Machida fc Hosokawal (l2013ll . which include magnetic 


ope that ancho rs the 


age ((Yen et al 


2914 . The 


helds, have successfully demonstrated the growth of protostars and disks, as well as the 
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deceleration of infall and the dissipation of envelopes in the late evolutionary stage. Their 
models also predict that, in the last phase of evolutionary stage, the g rowth of the centr al 


Keplerian disk stops while t he protostellar mass itself is still growing fjVisser et ah 


Machida & Hosokawa 

2013; 

Machida et ah 

2014) 


2009 


2014J) . Currently our observational samples do 


not show any clear evidence for such a saturation of the disk radius. Future large-scale 
ALMA surveys of a number of protostellar sources may reveal such a saturation, and a 
comprehensive observational picture of disk formation and growth around protostellar 
sources. 


6. Conclusions 

Using ALMA in Cycle 0, we observed a Class I protostar in Taurus star-forming region, 
TMC-IA in the 1.3-mm dust continuum, the ^^CO (J = 2 — 1), and C^®0 (J = 2 — 1) lines. 
The main results are summarized below. 

1. The 1.3-mm dust continuum shows a strong compact emission with a weak extension 
to the west. The gas mass of the dust emission is estimated from its flux density to be 
Afdust = 4.2 X 10“^ Mq assuming a opacity Ki.smm = 1-2 cm^g“\ Tc = 28 K and the 
gas/dust mass ratio of 100, which is ten times larger than the gas mass estimated from the 
C^®0 flux. 

2. The ^^CO line traces the molecular outflow with the axis perpendicular to the 
elongation of the 1.3-mm continuum emission. Velocity gradients are seen both along and 
across the outflow axis. The velocity is accelerated along the outflow axis (Hubble law). 
Across the outflow axis, the higher velocity component is located closer to the axis, which 
can be explained by an outflow driven by a parabolic wide-angle wind. 

3. The C^®0 emission shows the velocity gradient along the major axis of the disk 
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traced in the 1.3-mm continuum emission. This velocity gradient is due to rotation. The 
power-law indices of the radial prohle of the rotation velocity are estimated from httings 
to the major-axis PV diagram to be Pin = 0.54 for the inner/high-velocity component and 
Pout = 0.85 for the outer/low-velocity component. This indicates the existence of the inner 
Keplerian disk surrounded by the outer infalling envelope. 

4. In order to investigate the nature of the Keplerian disk, httings are performed to 
the channel maps at high velocities with a rotating disk models. In the model htting, 
the power-law index of the rotational velocity (n) is included as a free parameter. The 
reduced is minimized at n = 0.53, which is very consistent with pi^, and conhrms that the 
rotational velocity at high velocities can be explained by Keplerian rotation. The dynamical 
mass (not inclination corrected) derived from our best-ht disk model is M* sin^ i = 0.56 Mq. 

5. In addition to the velocity gradient along the major axis, there is another velocity 
gradient along the minor axis, which can be interpreted as infalling motions. In order to 
investigate the nature of the infall motion, the observed PV diagrams along both major and 
minor axes are compared with models consisting of an infalling envelope and a Keplerian 
disk. The Keplerian-disk size and the inclination angle are estimated to be 100 AU and 

i = 65° respectively, and the infall velocity is found to be 0.3 times as large as the free fall 
velocity yielded by the dynamical mass of the central protostar. Possible reasons for the 
infall motions to be slower than the free fall velocity are magnetic pressure and magnetic 
tension. The magnetic held strength yielding the slow infall velocity is estimated to be ~ 2 
niG. 

6 . Based on observations, including our series of ALMA observations, which show 
transitions from infalling envelopes to Keplerian disks, evolutionary processes of Keplerian 
disks around protostellar sources are examined. Those samples exhibit correlations in Tboi 
vs M*, Lfooz/-hsubmm VS M*, and M* vs i?kep- This may suggest that the central protostellar 
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masses and disk radii increase as protostellar evolution proceeds. 

This paper makes use of the following ALMA data: ADS/JAO.ALMA^2011.0.00210.S 
(P.I. N. Ohashi). ALMA is a partnership of ESO (representing its member states), NSF 
(USA) and NINS (Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan), in 
cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, 
AUI/NRAO and NAOJ. We thank all the ALMA staff making our observations successful. 
We also thank the anonymous referee, who gave us invaluable comments to improve the 
paper. 

Facilities : ALMA 


A. Analytic Consideration on our Method to Derive Rotational Profiles from 

Observed PV Diagrams 


In section IT^ we adopt two independent methods to derive the rotational profile of the 
disk around T MC-lA; one is a method using the observed PV diagram along the major axis 


introduced by 


Yen et ah 


(120131) . and the other a full three-dimensional htting of a rotating 


disk model to the observed C^®0 channel maps in the high-velocity (> 2.4 km s“^) regions. 
The former method yields the rotational profile of 2.4 (km s“^) x {R/67 and 

the latter 2.4 (km s“^) x {R/87 Thus, while the derived rotational power-law 

indices are consistent with each other and suggest Keplerian rotation, the rotational radius 
at a certain velocity derived from the former method is ~ 23% lower than that from the 
latter method. This implies that the former method underestimates the central stellar mass 
and the rotational angular momentum by ~ 23%. In this appendix, these uncertainties 
inherent in the method using the PV diagrams are analyzed. For simplicity, we consider a 
geometrically-thin, uniform disk with a pure rotational motion, and discuss the effects of 
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the finite velocity and spatial resolutions. The effect of the internal gas motion in the disk 
can be absorbed into the effect of the finite velocity resolution. We do not consider any 


radial motion in the disk here, since the effects of 
rotational profiles have already been discussed by 


he radial motions on the estimates of the 


Yen et ah 


(120131). 


The on-plane coordinates of the model disk are described with the 2D polar coordinates 
system {R,(j)), and the projected coordinates on the plane of the sky 2D cartesian {x,y), 
where the x-axis is defined as the major axis of the ellipse of the disk projected on the 
sky. The disk inclination angle from the plane of the sky (= i) is defined such that i = 0 
means face-on, and thus x = Rcoscj) and y = i?sincosh The radial profile of the rotation 
velocity of the disk is expressed as V^{R) = Vi{R/Ri)~'^, and the iso-velocity curve at a 
given line-of-sight velocity (= 14os) projected onto the plane of the sky is described as 
cos0sink This curve passes through the origin {x,y) = (0,0) and (i?o, 0) 
where Ro = Ri{Vi^os/Vi/ sinz)”^/'^. Ro gives the outermost point located in the major axis 
at a given 14os and curves at different Wos do not intersect with each other. The curves 
exhibit a symmetrical shape with respect to the major axis, and around the systemic velocity 
the curves show the well-known “butterfly” shapes. The curves at different 14os correspond 
to the velocity channel maps in the case of the infinite spatial and velocity resolutions, 
and the PV diagram along the major axis forms the rotation curve of Wos = h'^(-R) sinz 
exactly. In cases of finite spatial and velocity resolutions, the iso-velocity curve on the 
plane of the sky must have widths depending on the velocity resolution (= AD), where 
the maximum width is seen along the major axis. Such iso-velocity regions are convolved 
with the spatial beam, defined as the circular Gaussian with the FWHM beam size of 6 
(= exp[—41n2(a;^ -|- y‘^)/6‘^]). The spatial convolution skews the mean position along the 
major axis {= Xm) at a given Wos- 

In realistic cases of central disks, our method adopts only high-velocity (typically 
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> a few times km s“^), fast-rotating parts of the P-V diagrams, and the amount of the 
positional shifts within a typical velocity resolution (~ 0.1 km s“^) is much smaller than 
a typical spatial resolution. The iso-velocity curve (ie., Id^os = V 0 (i?) cos 0sini) can be 
re-written in terms of y as 


Ro cos i 


y 



2/h+i) 



2 


(Al) 


Thus, the width of the iso-velocity region along the y direction at a given x{< Ro) and Al/ 
can be written as 



(A 2 ) 


where (A1//1 /los)^ or higher order components are ignored. On the major axis {y = 0), the 
intensity distribution (= I) spatially convolved along the y axis is then described as; 



2 


(A3) 


le^ {z ), oc Ay exp -4 In 2 


where ^ = x/Ro is the normalized spatial coordinate along the major axis, and 

= (^/Ro/ cosi is the normalized beam width along the minor axis. In Eq. flA3p . Ay‘^ or 
higher order components are ignored because Ay is an inhnitesimal proportional to AV. 
The spatial convolution along x-axis (the major-axis direction) does not change the mean 
position mathematically. Then, the skewness of the mean position along the major axis due 
to the hnite resolution AV and 9 can be expressed as; 



(A4) 


where Zm is the normalized mean position. Zm can be derived as a function of 9n, which 
depends on Ro (or 1 /los) because 6 * is a constant. Fig. ITTI shows the dependence of Zm on 
l/9n in the case of n = 1/2 (Keplerian rotation). This hgure shows that large l/9n gives 
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Zm — 1, and thus the relevant radius can be derived accurately if the spatial resolution 
is high enough as compared with the emission size. On the other hand, small l/6*„ gives 
Zm — 0.760 < 1, which means that the correct radius cannot be derived from the mean 
positions. However, if the beam size is very large {l/9n < 1.0), the dependency of Zm on 
l/6n (i.e. Ro) is small, indicating Xm = const, x Rq. This means the power-law index 
can be derived correctly by using the mean positions, though the radii themselves are 
underestimated. These results also indicate that the estimates of the central stellar masses 
or angular momenta are affected at most <23%. 

In our case of TMC-IA, the size of the synthesized beam (FWHM~ 1") corresponds 
to l/9n — Ro/^'^ = i7o/200 AU with i = 65° and d = 140 pc. Fig. [T4l shows that the 
power-law index of the rotational velocity is expected to be estimated correctly at i? < 200 
AU for TMC-IA. In an intermediate situation, the degree of the underestimation of the 
radii depends on the relevant radius, thus the power-law indices must also be distorted. 
These considerations imply that effects of the spatial and velocity resolutions, as well as 
those of the radial velocities, should be taken into account to derive the rotational prohles 
from the observed PV diagrams. 


B. Comparison between the flattened and spherical envelope models 


In Sec. 15.21 we have constructed models consisting of a Keplerian disk and a flattened 
envelope (FE model) to interpret the observed gas motions, and have found that the 
envelope is infalling toward the disk with the infalling velocity slower than the free fall 
velocity yielded by the mass of the central star. The interpretation of gas motions are . 


however, dependent on the assumed geometry of the envelope (e.g. 


Tobin et ah 


2 ni2ah . 


and it is thus possible that the identihed infalling motion slower than the free fall could be 
due to the effect of the assumed envelope geometry. To qualitatively assess the effects of 
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the different geometries of the envelopes, we have also constructed models consisting of a 
spherical envelope and a Keplerian disk. Hereafter we call them SE models. 


The density structure of the spherical envelope in the SE models we adopt ed is a 


solution of the rotating and infalling envelope modeled by 


Terebev et ah 


( 1984 ). with 


the temperature prohle of T(r) = 38.0 K (r/100 the mass infalling rate of 

M = 10“^ Mq/jt (M*/0.5 the central stellar mass M* = 0.68Mq, and the 

outermost radius R = 900 AU. To mimic the cavity created by the associated outflow, 
artihcial bipolar cavities in a conical shape with the opening angle of 20° are added in the 
spherical envelope, and within R < i?kep = 100 AU the same disk as that of EE models is 
located. To simplify the comparison of the observed gas kinematics between EE and SE 
models, the rotational velocity in the spherical envelope is assumed to be only the function 
of the distance from the rotational axis, with the same radial prohle as that of the battened 
envelope models, and the infalling motion is assumed to be isotropic. As the case of EE 
models, we adopt two types of infalls in SE models; a reduced infall K- = 0.3Ufr, and 
free fall (K = Uff), where the free fall velocity is derived from the central stellar mass of 
0.68 Mq. The former and latter models are named as SE model A and D, respectively. The 
corresponding EE models are EE model A (our best ht model) and D. 


Figure [T5h and b compare the PV diagrams along the major and minor axes of EE 
model A (red contours) and SE model A (grey scale), and Fig. [TR ; and d compare EE 
model D (red contours) and SE model D (grey scale), respectively. These hgures show some 
slight differences between SE and EE models as described below. As shown in Fig. [TM . the 
minor-axis PV diagram of SE model D shows two peaks at U ~ 5 km s“^ and ~ 8 km s“^, 
while two velocity peaks in EE model D are ~ 1 km s“^ faster than them. In addition, 
EE model A/D shows only blueshifted emission at positions ^ —1" in the minor-axis PV 
diagram and only redshifted emission at 1", whereas SE model A/D shows blue- and 
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redshifted emission on both positive and negative sides. In spite of these slight differences, 
FE model A and SE model A, as well as EE model D and SE model D, show overall similar 
PV diagrams along both major and minor axes and we, therefore, conclude that SE models 
cannot explain the observations better than FE models. It is important to note that when 
SE model A and D are compared with the observations, it is found that SE model A can 
explain the observations better than SE model D, suggesting that the infall velocity should 
be significantly slower than the free fall velocity even in the case with a spherical envelope. 
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J2000 Right Ascension 


Fig. 1.— Continuum emission map of TMC-IA observed with ALMA. Contour levels are 
—3, 3, 6,12, 24,... X(T where la corresponds to 0.96 mJy beam“L A cross shows the position of 
the continuum emission peak. A hlled ellipse at the bottom left corner denotes the ALMA 
synthesized beam; I'.'Ol x 0'.'87, P.A. = +0.87°. The deconvolved size of the continuum 
emission is Ot'bO ± 0703 x 0737 ± Ot'OS, P.A = +73° ± 16° and the elongation direction (73°) 
is shown with a white dashed line. Blue and red arrows show the direction of the molecular 
outflow from TMC-lA observed in ^^CO (J = 2 — 1) line (Fig. [2]). 
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J2000 Right Ascension 


Fig. 2.— Moment 0 {white contours) and 1 {color) maps of the (J = 2 — 1) emission in 
TMC-IA. Contour levels of the moment 0 map are —3, 3, 6,12, 24,... xa where la corresponds 
to 0.034 Jybeam“^km s“^. A Cross shows the position of the central protostar (continuum 
emission peak). A hlled ellipse at the bottom-left corner denotes the ALMA synthesized 
beam; 1"02 x 0"90, P.A. = —178°. ^^CO clearly traces the molecular outflow from TMC-lA 
and the axis of the outflow (~ —17°) is shown with a dashed line. 
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velocity (km/s) 


Fig. 3.— Position-Velocity diagrams of the (J = 2 — 1) emission in TMC-lA along the 

outflow axes (PA = —17°). Contour levels are 6a spacing from 3a {la = 20 mJybeam”^). 
Central vertical dashed lines show the systemic velocity (14ys = 6.4 km s“^) and central 
horizontal dashed lines show the central position. 
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J2000 Right Ascension 


Fig. 4.— Moment 0 {white contours) and 1 {color) maps of the (J = 2 — 1) emission in 
TMC-IA. Contour levels of the moment 0 map are —3, 3, 6,12, 24,... x a where la corresponds 
to 8.1 mJybeam“^km s“^. A cross shows the position of the central protostar (continuum 
emission peak). A hlled ellipse at the bottom-left corner denotes the ALMA synthesized 
beam; 1'.'06 x (/.'QO, P.A. = —176°. Blue and red arrows show the direction of the molecular 
outflow from TMC-lA observed in (J = 2 — 1) line (Fig. [2]). 
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Fig. 5.— Channel maps of the C^®0 (J = 2 — 1) emission in TMC-lA. Contour levels are 
—3, 3, 6,12, 24,... X a where la corresponds to 7.1 mJybeam”^. Crosses show the position of 
the central protostar (continuum emission peak). A hlled ellipse in the top left panel denotes 
the ALMA synthesized beam; 1"06 x 0'.'90, P.A. = —176°. LSR velocities are shown at the 
top-left corner of each panel and the systemic velocity is 6.4 km s“^. 
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Fig. 6.— Maps of blueshifted (blue contours) and redshifted (red contours) emission of (a) 
high-velocity (|A1/| > 2.0 km s“^) and (b) low-velocity (|A1/| < 1.0 km s“^) components 
of (J = 2 — 1) emission in TMC-lA as observed with ALMA. Contour levels are 

—2,2,4,8,16,32 x a where la corresponds to 2.9 mJybeam”^ and 4.1 mJybeam”^ for (a) 
and (b), respectively. Crosses in both panels show the position of the central protostar 
(continuum emission peak). Filled ellipses at the bottom left corner in both panels are the 
same as Fig. 01 The dashed line passing peaks of blueshifted and redshifted maps in each 
panel shows the direction of the velocity gradient. 
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offset [arcsec] 

Fig. 7.— 2D mean positions of the (J = 2 — 1) emission in TMC-lA. Each point 

indicates a 2D mean position at each channel weighted with the intensity and 6a cutoff 
{a is the rms noise level of emission). Error bars are not plotted but are smaller 

than the point size. Dashed lines show the major and the minor axis of the disk (73°, —17° 
respectively) and the intersection is the position of the central protostar (continuum emission 
peak). The point size becomes smaller as the value of |17 — 14ys| becomes higher and green 
points correspond to a velocity range of |A1/| > 2.4 km s“^. 
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Fig. 8.— Position-Velocity diagrams of the (J = 2 — 1) emission in TMC-lA along 

(a) the major and (b) the minor axes (the major axis corresponds to the white dashed line 
in Fig. [H PA = 73°). These PV diagrams have the same angular and velocity resolutions 
as those of the channel maps shown in Fig. [5l Contour levels are 6a spacing from 3a 
{la = 7.1 mJybeam”^). Central vertical dashed lines show the systemic velocity (Vys = 
6.4 km s“^) and central horizontal dashed lines show the central position. Points with error 
bars in the panel (a) are mean positions derived along the position (vertical) direction at 
each velocity. Curves in the panel (b) show escape velocities along the outflow-axis direction 
for i = 70° and 40°. 
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Fig. 9.— Mean positions of the PV diagram along the major axis plotted on logi?—logl/ 
plane. The ordinate is not de-projected. Blue and red points show blueshifted and redshifted 
components, respectively. Dashed lines show the best-£t lines with a double power law. In 
addition to the inner and outer power, the “break” radius and velocity are included as free 
parameters. The best £t parameters are Rb = 67 AU, Vj, = 2.4 km s“\ pin = 0.54 ± 0.14, 
Pout = 0.85 ± 0.04. The error for {Rb, H) is ~ 1%. 
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Fig. 10.— Channel maps, nsed in a fitting with the rotating disk models of (a) the best-fit 
disk model (red contonrs) snperposed on that of the observations of (J = 2 — 1) emission 
and (b) the residnal obtained by snbtracting the best-fit disk model from the observations. 
The systemic velocity is 14sr = 6.4 km s“^. Contour levels are —3, 3, 6,12, 24 ,... x a where 
Icr corresponds to 7.1 mJybeam”^ for the panel (a) and —3, —2, 2, 3,4... x a for the panel 
(b). The spatial scale is different from Fig. [5l 
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Fig. 11.— Distribution of the reduced terms of the inclination angle of the disk. All 


points have the same parameter set as the best-fit disk model except for the inclination 
angle. Larger inclinations (45° < i < 75°) give small reduced (that is, good) models to 


explain the observations. 
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Fig. 12.— Comparisons of the observations and models with infall motions. The left column 
shows the mean position of the PV diagram along the major axis in log R — log V diagrams. 
Red points indicate models while black points indicate the observation, which is the same 
plot as Fig. [9] except for the color. The middle and right columns show PV diagrams along 
the major and minor axis respectively. The observations are in black contours and grayscales 
while models are in red contours. Contour levels of the PV diagrams are the same as Fig. 
[HI The same row includes the same model and the parameters of each model are indicated 
in the left panel of each model. FE model A explains the observations most reasonably. 
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Fig. 12.— (Continued) 
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Fig. 13.— Correlation between (a) bolometric temperatures and central protostellar masses, 
(b) Lboi/Lsnhram and Central protostellar masses, and (c) central stellar masses and sizes of 
Keplerian disks around protostars. The data values are listed in Table [3l Red filled circles 
corresponds to the protostars for which the transitions from infall motions to Keplerian 
rotations are identihed well. 
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Fig. 14.— The dependency of the normalized mean position Zm = Xm/Ro on the normalized 
beam width l/9n = RoCosi/6 in the case with v = 1/2 (Keplerian rotation). 
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Fig. 15.— Comparison FE model A/D in Sec. 15.21 (red contours) and SE model A/D in 
Appendix [B] (black contours and grayscales). The right and left columns show PV diagrams 
along the major and minor axes, respectively. The contour levels are the same as Fig. [Hland 
Fig. m 
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Table 1. Summary of the observational parameters 


Interferometer and date 


CO J = 2 — 1 J = 2 — 1 continuum 

ALMA (extended), 2012.Nov.6 


Target 

Coordinate center 

Frequency 
Primary beam 
Projected baseline length 
Synthesized beam (P.A.) 

Velocity resolution 
Noise level (no emission) 
Noise level (detected channel) 
Passband calibrator 
Flux calibrator 
Gain calibrator 


TMC-IA 

R.A. (J2000)=4*'39“'35=010 
Dec. (J2000)=25°41'45''500 


230.5380 GHz 

219.5604 GHz 

225.4336 GHz 

27'.'3 

28'.'6 

15.9 - 285.7 m 

2779 

'02 X 0"90 (-178°) 

1706 X 0790 (-176°) 

1701 X 0787 (+0.87‘ 

0.16 kms”^ 

0.17 kms-3 

923 MHz 

5.6 mjy beam-3 

4.6 mjybeam—3 

— 

20 mjybeam—3 

7.1 mjybeam—3 

J0522-364 

Callisto 

J0510+180 

0.96 mjy beam-3 


Table 2. Fixed and free parameters of the model htting 


fixed 

i 

65° 

Riu 

0.1 AU 

-Rout 

200 AU 



free 

best 

0.56t°;°s Mq 

M 200 *' 

3.6ti o 10“^ ^0 

p 

l.46t?;3 

Tioo 

38t® K 

Q -Rcent 

-0.02l3;°l 166t«° AU 0.53t° 


= M* sin^ i, not inclination corrected. 

'^The mass within Rout = 200 AU. The mass of the Keplerian disk (iikep = 100 AU) is calculated to be 
2.5 X 10-3 Mq. 










Table 3. Parameters of Keplerian disks around protostars 


Source 

Lbol (Lq) 

Tboi (K) 

^bol! -^submm 

A^kep (AU)^ 

M* (Mq^ 

Afdisk (Af©)'* 

Class 

transition?® 

References^ 

NGC1333 IRAS4A2 

1.9*5 

51*5 

28*5 

310 

0.08 


0 

no 

1 

VLA1623A 

1.1 

10 

83 

150*' 

0.22 


0 

yes 

2,13 

L1527 IRS 

1.97 

44 

45 

54 

0.30 

0.0028-0.013 

0 

yes 

3,12,16 

R CrA IRS7B 

4.6 

89 


50 

1.7 

0.024 

I 

no 

4 

L1551 NE 

4.2 

91 

111 

300 

0.8 

0.026 

I 

yes 

5,14,17 

L1551 IRS 5 

22.1 

94 

107 

64 

0.5 

0.07 

I 

yes 

6,12,16 

TMCl 

0.9 

101 

114 

100* 

0.54 

0.025-0.06 

I 

no 

7,12,16 

TMC-IA 

2.7 

118 

143 

100 

0.68 

2.5 X 10-3 

I 

yes 

8,12,16 

TMRl 

2.6 

140 

734 

50 

0.7 

0.01-0.015 

I 

no 

7,16 

L1489 IRS 

3.7 

238 

259 

700 

1.6 

3-7x10-3 

I 

yes 

9,15,16 

L1536 

0.4 

270 


80 

0.4 

0.07-0.024 

I 

no 

7,12 

Elias 29 

14.1 

299 

4215 

200 

2.5 

< 7 X10-3 

I 

no 

10,12,16 

IRS 43 

6.0 

310 


140 

1.0 

8.1 X 10-3 

I 

no 

11,15 

IRS 63 

1.0 

327 

33 

100 

0.37 

0.055 

I 

no 

10,12,16 


^^submm is defined as the luminosity measured at wavelengths longer than 350 /rm. 
Outer radius of the Keplerian disk. 


‘^Central protostellar mass. 

"^Mass of the Keplerian disk. 

®Frag showing whether or not sources show transitions from infall motions to Keplerian rotations, “yes” and “no” indicate sources with and 
without transitions, respectively. 






^References: fli Choi et al. 

2010 

; (2) 

Murillo et al. 2013 


Ohashi et al. 

2C 

)14 


Lindbere et al. 2 

014 

; (5) 

Takakuwa et al. 

2014 

; (6) 

Chou et al. 

2014 

; (7) 

Harsono et al. 

2014 

; (8) this work: (9) 

Yen et al. 

2014 

; (10) 

Lommen et al. 

2008 

; (11) 

(17) 

Eisnei 

2012 

; (12) 

Kristensen et al. 

201S 

; (13) 

Murillo & Lai 

2013 

(14) 

FVoebrich 2005: fl5i Joreensen et a 

. 2009; 

(16) 

Green et al. 

2013 

Proebrichl 

2005 


SBoth Lbol Q-nd -Lsub assumed to be a half of those derived toward IRAS 4A binary system. Thoi is also assumed to be the same as that 
of the binary system. 


Murillo et al 


I 2013) concluded that the Keplerian radius of VLA1623A is 150 AU though their UV-space PV diagram shows a turn over 
point of the rotational velocity profile at 50 AU. They argued that Rkep t)® larger than 50 AU because of optical depth and absorption. 

^The analysis using logR — logU diagram gives us a breaking radius of ~ 90 AU. The results here is, however, not affected significantly even 
if 90 AU is adopted instead of 100 AU. 
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